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ABSTRACT 



Aims. We analyze OGLE-2007-BLG-050, a high magnification microlensing event (A ~ 432) whose peak occurred on 2 May, 2007, with 
pronounced finite-source and parallax effects. We compute planet detection efficiencies for this event in order to determine its sensitivity to the 
presence of planets around the lens star. 

Methods. Both finite-source and parallax effects permit a measurement of the angular Einstein radius 9e = 0.48 ± 0.01 mas and the parallax 
tte = 0. 12 ± 0.03, leading to an estimate of the lens mass M = 0.50 ± 0.14Mo and its distance to the observer = 5.5 ± 0.4 kpc. This is only 
the second determination of a reasonably precise (< 30%) mass estimate for an isolated unseen object, using any method. This allows us to 
calculate the planetary detection efficiency in physical units (r±,m,,), where r± is the projected planet-star separation and m,, is the planet mass. 
Results. When computing planet detection efficiency, we did not find any planetary signature and our detection efficiency results reveal signif- 
icant sensitivity to Neptune-mass planets, and to a lesser extent Earth-mass planets in some configurations. Indeed, Jupiter and Neptune-mass 
planets are excluded with a high confidence for a large projected separation range between the planet and the lens star, respectively [0.6 - 10] 
and [1.4 - 4] AU, and Earth-mass planets are excluded with a 10% confidence in the lensing zone, i.e. [1.8 - 3.1] AU. 



Key words, extrasolar planets - gravitational microlensing 
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1. Introduction 



Over the last decade, microlensing events have been intensively followed in order to detect extrasolar planets around lens stars 
and to measure their abundance in our Galaxy. This is one of the few planet-detection techniques that is sensitive to very low 
mass planets, and microlensing discoveries have twice the record for the lowest mass planet to orbit a star other than a stellar 
remnant (Beaulieu et al. 2006; Bennett et al. 2008). During a microlensing event, i.e. when a background source passes close to 
the line of sight to a foreground lens star, the observed source flux is magnified by the gravitational field of the lens. The presence 
of a companion around the lens star introduces two kinds of caustics into the magnification pattern : one or two "planetary 
caustics" associated with the planet and a "central caustic" close to the primary lens projected on the sourc e plane. When the 



source crosses or approach es one of these features, deviations appear from a single point-lens light curve ao & Paczvnskil 
199lHGould&Loeb(ll992h . 



1.1. Central caustic and detection efficiency 

Signific ant efibrt has been expended on the observation and modeling of hig h magnification events because they probe the central 

Any planets in the system are highly likely to affect 



aignmc ant errort has been expended on the observation and modelmg or hig 
caustic dOriest & Safizadehlll998l:lRhie et allboOOl: iRattenburv et al.li2002h . 



the central caustic, resulting in potentially high sensitivity to the presence of even low-mass planets. 

Indeed, a major advantage of the central caustic is that it is possible to predict in advance when the source passes close to 
the line of sight of the lens and so when there is the greatest chance of detecting planets. Thus observations can be intensified, 
further improving the sensitivity to planetary-induced anomalies in the lightcurve. 

In these specific cases, for which the impact parameter can be very small, finite -source effects might strongly affect and 
diminish a possible planetary signal (e.g.. iDong et al. 2009b , Bennett et Rhie |[T996h . In the absence of any deviation from a 
finite-source single point-lens model, one can still compute the planet detection efliciency in order to derive upper limits on the 
probability that the lens harbors a planet (Gaudi & Sackett 2000). It also allows to combine statisticall y the detection efficiencies 
computed from observed events to estimate the frequency of planetary companions to the lens dGaudi et alJi2002l) . 

The extremely high magnification microlensing event OGLE-2007-BLG-050 was well followed and is a good candidate for 
analyzing the sensitivity of such an event with pronounced finite-source effects to the pre sence of a planeta ry companion. In this 
study, we compute the planetary detection efficiency for this event, following the iGaudi & Sacketi tood) method. To perform 
the calculations of binary light curves, we use the binary-lens finite-source algorithm developed by lDong et al. I (l2006l) and the 
formalism of .Yoo et al.. (.2004a) for the single-lens finite-source effects. 



1.2. Mass and distance estimates of the lens star 

OGLE-2007-BLE-050 is also one of the rare events that can potentially be completely solved by measuring both the microlens 



Einstein angular radius (9e and the microlens parallax tte- Indeed, after the first microlenses were detected (iAlcock et al.llI993 
[Udalski et al..,1993,) . several authors showed that the microlens Einstein angular radius 6e, 



could be measured from deviations relative to the standard point-lens ([Paczvnski Il98d) lightcurve, due to finite-source effects 
( Gould et al.ll994 ; Nemiroff" & Wickramasinghell994 ; Witt & Maoll994 ). The measured parameter associated with these effects 
is p., corresponding to the angular size of the source 0, in units of Oe - The measur ement of 6e constrains the physical properties 
of the lens and so leads to the first part of a full solution for an event (lGouldll2000h . 



6e = -'jKMnri 



AG 

c^AU 



8 mas M„ 



where M is the lens mass and n^ei is the lens-source relative parallax. For most events, the only measured parameter that depends 
on the mass M is the Einstein timescale, fe, which is a degenerate combination of the lens mass M, the lens-source relative 
parallax ;7r,ei and the proper motion /irei. It can be expressed as : 



Gouldl (119921) showed that if one measures both 6e and the microlens parallax, tte, which is derived from the distortion of the 



microlens light curve induced by the accelerated motion of the Earth, one can determine 



TTe = 




* Royal Society University Research Fellow 
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and so determine the lens mass and lens-source relative parallax as well, 

M - ; JTrel = TTeOe- 

KRe 

After thousands of single-lens microlensing events discovered to date, measurements of both Op, and tte still remain a chal- 
lenge. The mic rolens parallax tte has been measured for more than twenty single lenses (lAlcock et al ]|i995 [the first parallax 
measurement!. IPoindexter et alJ 2005 a nd ref erences therein), while the angular Einstein r adius 0r has been measured for only 
few cases of single lenses jAlcock et all 19971 200 1 ; ISmith et alj|2003bl ; IYoo et alj|2004al ; Ijiang et alj|2004l ; ICassan et alj|2006l 
Gould et al. 200^. 



However, reliable m ass estimates for isolated stars have been determined with microlensing only twice. lAlcock et al.l (1200 lb 
and lGould etaP (l2009l) each measured both Oe and tte respectively for MACHO LMC-5 and OGLE 2007-BLG-224. For MACHO 
LMC-5, good measurements of TTrei and //,ei were obtained with the original photometric data and additional high resolution 
photometry of the lens (HST observations). Only for OGLE 2007-BLG-224 has there been a reliable mass estimate derived using 
only ground-based photometric data. 

All other good microlens stellar mass measurement s to date have been obtained for binary (or pl anetary) lens events : EROS 
BLG-2000-5 (lAnetal.ll2002h. OGLE 20 06-BLG-109 (iGaudi et al. 20 08b, OGLE 20 07-BLG-071 dPong et alJ l2 009al). OGLE 
2003- BLG-267 Jjaroszynski et al.ll2005l) . OGLE 2002-BLG-069 (.Kub"as et al.ll2005i) and OGLE 2003-BLG-235 JBond et af 



20041). 



1.3. Detection efficiency in physicai units 

Here, we present ground based photometric data of the event OGLE 2007-BLG-050 which we use, for the first time, to constrain 
both the presence of planets and the mass of the lens. 

This is also the first event for which parallax and xallarap (source orbital motion) are analyzed simultaneously. However, we 
find that the apparent xallarap signal is probably due to minor remaining systematic effects in the photometry. 

Access to the physical properties of the lens allows us to compute the planetary detection efficiency in physical units (r^, nip), 
where rj_ is the projected separation in AU between the planet and the lens and nip is the planet mass in Earth mass units. 

OGLE-2007-BLG-050 had a high sensitivity to planetary companions of the lens, with a substantial efficiency to Neptune- 
mass planets and even Earth-mass planets. 



2. Observational data 



The microlensing event OGLE-2007-BLG-050 was identified by the OGLE III Early Warning System (EWS ; lUdalskill2003h 
(a = 17/!58ml9.39,s, 6 = -28°38'59" (J2000.0) and / = +1.67", b = -2.25°) on 2 Mar 2007, from observations carried 
out with the 1.3 m Warsaw Telescope at the Las Campanas Observatory (Chile). The peak of the event occured on HJD' = 
HJD - 2, 450, 000 = 4221 .904 (2007 May 2 at 9:36 UT). 



The event was monitored over the peak by the Microlensing Follow Up Network (//FUN : lYoo et al ]|2004a') from Chile (1.3m 



SMARTS telescope at the Cerro Tololo InterAmerican Observatory), South Africa (0.35 m telescope at Bronberg observatory), 
Arizona (2.4 m telescope at MDM observatory, 1.0 m telescope at the Mt Lemmon Observatory), New Zealand (0.40 m and 
0.35 m telescopes at Auckland Observatory and Farm Cove observatory respectively) and on the wings from the Vintage Lane 
(Marlborough, New Zealand), Wise (Mitzpe Ramon, Israel) and Palomar 60-in (Mt Palomar California, USA) observatories. 
However, the last three were not included in the final analysis because they do not significantly improve the constraints on 
planetary companions. Data from all the three sites are consistent with single-lens model. 

It was also monitored by Microlensing Observations in Astrophysics (MOA) with the 1.8 MOA-II telesco pe at Mt John 
University Observatory (New Zealand), and Probing Lensing Anomalies Network (PLANET ; Albrow et alJ 1998h from 5 differ- 



ent telescopes : the Danish 1 .54 m at ESO La Silla (Chile), the Canopus 1 m at Hobart (Tasmania), the Elizabeth 1 m at the South 
African Astronomical Observatory (SAAO) at Sutherland, the Rockefeller 1.5 m of the Boyden Observatory at Bloemfontein 
(South Africa) and the 60 cm of Perth Observatory (Australia). The RoboNet collaboration also followed the event with their 
three 2m robotic telescopes : the Faulkes Telescopes North (FTN) and South (FTS) in Hawaii and Australia (Siding Springs 
Observatory) respectively, and the Liverpool Telescope (LT) on La Palma (Canary Islands). 

In this analysis, we use 601 OGLE data points in / band, 104 //FUN data points in / band, 77 //FUN data points close to R 
band, 121 PLANET data points in / band, 55 RoboNet data points in R band and 239 MOA-Red data points (wide band covering 
R and / bands). 



3. Event modelling 

OGLE-2007-BLG-050 is a very high magnification event (A ^ 432) due to its small impact parameter uq. Because they are quite 
obvious on the observed light curve, finite-source effects must be incorporated in the modeling. Moreover, the long timescale of 
the event implies that parallax effects are likely to be detectable. 
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3.1. Finite-source effects 

When observing a microlensing event, the resulting flux for each observatory-filter / can be expressed as, 

Fi(t) = F,jA[u(t)] + Fhj, 

where , is the flux of the unmagnified source, Fi, , is the background flux and u(t) is the source-lens projected separation in the 

lens plane. 

When the so urce can be approximated as a point, the magnification of a single-lens event is given by dEinstein Il936l : 
Paczvhskilll986h 

A(m) = — — 

u Vm^ + 4 

However, in our case the source cannot be consider ed as a point (u < Pt) and the variation in brightness of the source star across 
its disk must be considered using the formalism of Yoo et al. ('2004a'). When limb-darkeni ng of the source profile are neglected 
(uniform source), the magnification can be expressed as (.Gould. 1994a; .Witt & Mao 1 994; i Yoo et al.ii2004a) . 



A„„/(m/p«) - A(m)Bo(m/p»), 



Bo(z)= -zE(k,z) 
n 



where E is the elliptic integral of the second kind and k - min(z^\ 1). Separating the u and z = m/p» parameters allows fast 
computation of extended-source effects. 

To include the limb-darkening, we parameterize the source brightness S by. 



S{0) 



= 1 -r 



1 - -(1 - cos 61) 
2 



where 6 is the angle betwee n the normal to th e stellar surface and the line of sight. The new magnification is then expressed by 
adding the Bi(z) function of lYoo et al. I (l2004ah related to the linear limb-darkening law 

Au{u/p*)^A(u)[Bo(z)-rBi(z)]. 

The limb-darkening coeffic i ents F hav e been taken equal to 0.49 for the I filter and 0.60 for the R filter, which are results from a 
single-lens fit. From lClaret ( 2000t) and Afonso et al. (l2000l) models, considering a subgiant similar to our source (log g = 4, T = 
5250K), we find 0.44 and 0.53, respectively for I and R filters. These values are close to those of our model and lead to essentially 
the same parameter values as shown in TABLE 1 . 

In FiglT] we present the OGLE-2007-BLG-050 light curve modeled with extended-source effects (black curve) and without 
these effects (red curve). Finite-source effects are clearly noticeable by a characteristic flattening and broadening of the light 
curve at the peak. 

For each data set, the errors were rescaled to make per degree of freedom for the best-fit extended-source point-lens 
(ESPL) model close to unity. We then eliminated the largest outlier and repeated the process until there were no 3 cr outliers. 
None of the outliers constitute systematic deviations that could be potentially due to planets. 



3.2. Source properties from color-magnitude diagram and measurement of 9e 

To determine the dereddened color and magnitude of the microlensed source, we put the best fit color and magnitude of the 
source on an (/, V - 1), calibrated color magnitude diagram (CMD) (cf. Fig|2]). We use calibrated OGLE-III data. The magnitude 
and color of the target are / = 18.21 + 0.03 and {V - I) - 2.32 + 0.01. The mean position of the red clump is represented by an 
open circle at (/, V - I)rc - (15.95, 2.37), with an error of 0.05 for both quantities. The shift in position of our target relative to 
the red clump is then Al = 2.26 ± 0.05 and A(V - /) = -0.05 + 0.05. 



For the absolute clump magnitude, we adopt the Hipparcos clump magnitude Mj ^c = -0.23 + 0.03 dStanek & Garvanich 



(119981) ). The mean Hipparcos clump color of ( V-/)o,rc = 1 05+0.05 is adopted (Jennifer Johnson, 2008, privat e communication). 



Assum ing that the source is situated in the bulge and a Galactic center distance of 8kpc, /Jcc - 14.52 + 0.10 ( iEinsenhauer et al 



20051) 



The magnitude of the clump is given by /o,fic = A^/,«c + /"cc = 14.29 ± 0.10. We derive (/, V - /)o,rc - (14.29, 1.05) + 
(0.10,0.05). Hence, the dereddened source color and magnitude are given by : (/, V - /)o = A(/, V - I) + (I,V - I)o,rc - 
(16.55, 1.00)+ (0.12,0.08). 



From {V - /)o, we derive {V - K)o using the lBessel & Bretti ( 119881) diagram for giants, supergiants and dwarfs : {V - K)o - 



2.31 ± 0.13. The measured values of /q and (V - I)q then lead to Kq = 15.24 + 0.09. 

For completeness, we also derive an extinction estimate [Ai,E{V - I)] = (1.66, 1.32), which leads to an estimate Ryi - 
Av/E{V- I) ^2.02. 
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Fig. 1. Top: Light curve of OGLE-2007-BLG-050 near its peak on 2007 May \. Middle: Zoom onto the peak showing the finite- 
source effects. Bottom : Magnitude residuals. They correspond to the real residuals and are not exactly equal to the diff'erence 
between data and model of the light curve shown above, because the model is given in I band and the R band data points have 
been linearly converted into the I OGLE system. We show the model with finite source and parallax effects. As a comparison, a 
model without finite source effects is shown in red. 



The color determines the re lation between dereddened source flux and angular source radius. We use the following expression 
given bv lKervella et al.l (120041) for giants between AO and KO : 



log26». = 0.5170 - Q.lKii + 0.0755(y- K)o, 



giving 0^ = 2.20 + 0.06 nas. 

With the angular size of the source given by the extended source point lens (ESPL) fit, - 0.00458 + 0.00003, we derive the 
angular Einstein radius 9e'- 0^ - 6^1 p» - 0.48 + 0.01 mas, where the error is determined by : (crg^/OE)^ = + (o-pjpt)^. 

This first fit takes into account finite source effects only. The values of p, and 6e will not change significantly when adding new 
effects (see parallax effects later) but the induced modifications will be included in the final results. 

Then, combined with the fitted timescale of the event Ie = 66.9 ± 0.6 days, gives the geocentric relative lens-source proper 
motion : fi - Oe/Ie = 2.63 ± 0.08 mas/yr, with the same method for calculating the error. 



6 



V. Batista et al: OGLE-2007-BLG-050 : Planetary detection efficiency 



15 



1 \ r 



16 



17 



18 



19 



,»*„i-.C,j,.5j. 



20 



1.5 



2 

(V-l) 



2.5 



Fig. 2. Calibrated color-magnitude diagram of the field around OGLE-2007-BLG-050. The clump centroid is shown by an empty 
open circle, while the OGLE-III / and V - 1 measurements of the source are shown by an open circle surrounding 1 cr error bars. 



3.3. Parallax effects 

3.3.1 . Orbital parallax effects 

The source-lens projected separation in the lens plane, u(t) of Eq l3.1l can be expressed as a combination of two components, T(f) 
and y6(f), its projections along the direction of lens-source motion and perpendicular to it, respectively : 

u(t) = ^THt)+pHt). 

If the motion of the source, lens and observer can all considered rectilinear, the two components of u{f) are given by, 

t - to 

T(t) = — - ■ m = UQ. 

tE 

In the case of a simple point-source point-lens model, only five parameters are fitted : the source flux Fj, the blending flux Fb 
(both duplicated if more than one observatory), the time of the closest approach fo, the impact parameter mq and the timescale of 
the event fE- 

However, for long events, like OGLE-2007-BLG-050 (where Ie > yr/27r), the motion of the Earth cannot be approximated 
as rectilinear and generates asymmetries in the light curve. Parallax effects then have to be taken into account. To introduce these 
effects, we use the geocentric formalism (An et al. 2002 and Gould 2004) which ensures that the three standard microlensing 
parameters (fo, Ie, uq) are nearly the same as for the no-parallax fit. Now two more parameters are fitted. These are the two 
components of the parallax vector, tte, whose magnitude gives the projected Einstein radius, Te = AU/tte and whose direction is 
that of lens-source relative motion. 

The parallax effects imply additional terms in the Eq |3.3.1| 

T(f) = — + drit) ■ Pit) = Mo + m) 

tE 

where 

{6T{t),6/3{t)) = tte^s = (;rE.As,;rE x \s) 
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Fig. 3. Likelihood contours as a function of the parallax vector ^Te (1, 2, 3, 4 o"). The best fit is tie - (0.099, -0.072). There is a 
hard 3 cr lower limit tte > 0.086 which implies M < 0.67 Mq and a 3 cr upper limit tte < 0.23 which implies M > 0.25Mq. 



and As is the apparent position of the Sun relative to what it would have been assuming a rectilinear motion of the Earth. 

The Extended-Source Point-Lens (ESPL) fit yields a determination of the components {jie,n7'^e,e) of the parallax vector n^. 
projected on the sky in North and East celestial coordinates. This is done by mapping a grid over the tie plane and searching for 
the minimum of (cf. Fig|3]l. In addition to the best ESPL fit presented in Section [34l this grid search was done to probe the 
likelihood contours as a function of jte, holding each trial parameter pair tie - (nE,N, ^e,e) fixed while allowing all remaining 
parameters to vary. The best fit is t^e - (0.099, -0.072). There is a hard 3 cr lower limit tte > 0.086 and a 3 cr upper limit 
jte < 0.23. The error of tte is calculated from the 1 cr contour : ke - 0.12 + 0.03. The likelihood contours in the tte plane 
are slightly elonga ted along the North-South axis. This tendency, which is weak here due to the long timescale, is explained in 
Gould et aU (Il994{) by the fact that for short events the Earth's acceleration vector is nearly constant during the event. 

The Fig|4]shows the modeling improvement when we include the orbital parallax eff'ects in the fit. These plots only show the 
OGLE and MOA re siduals because thes e data mostly constrain the parallax since they cover a long time range. 

As discussed by lSmith et al.l (l2003al) . there is a mq <-> -mo degeneracy. For a low magnification event with \uo\ ~ 1, the uq > 
and UQ < solutions will behave dififerently, but for a high magnification event with |mo| « 1 like OGLE-2007-BLG-050, the 
Mo <-> -Mo transformation can be considered as a symmetry and there is no possibility to distinguish one solution from orbital 
motion alone. In principle, these can be distinguished from so-called "terrestrial parallax" effects caused by the different positions 
of the telescopes on the surface of the Earth. 



3.3.2. Terrestrial parallax effects 

We investigate terrestrial parallax in order to check if it is consistent with the vector parallax determined from orbital parallax 
eff^ects and to distinguish the mq > and mq < solutions. The resulting of the orbital-nterrestrial parallax model does not 
show any improvement and is actually worse than orbital parallax alone (Ax^ = 4, x\rbitai parallax ~ 1^60.5). The most likely 
explanation for this discrepancy is that the much stronger (Ax^ - "^^^^X without parallax ~ 1995.4) signal from orbital eff'ects reflects 
the true parallax and the small terrestrial parallax "signal" is actually just due to low-level systematic errors. 



3.3.3. Xallarap effects 



We also consider the possibility that the orbital parallax signal is actually due to xallarap (orbital motion of the source) rather 
than to real parallax. Of course an orbital motion of the source, in case of a binary orbit that fortuitously mimics that of the Earth, 
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Fig. 4. OGLE (stars) and MOA (hexagons) residuals (magnitude) for models with (upper panel) and without (lower panel) 
parallax effects. The residuals have been binned for clarity. 



can reproduce the same light curve as the orbital parallax effects but here we are looking for orbital motion that is inconsistent 
with the Earth-motion explanation. 

We therefore search for xallarap solutions (orbital motion of the source) by introducing 5 new parameters in the model related 
to the orbital motion of the source : P the period of the source's orbit, and ^e,e the xallarap vector which is analogous to the 

vector, and ai and 62, the phase and inclination of the binary orbit which function as analogs of the celestial coordinates of 
the source in case of parallax. The rather long timescale does not justify removing parallax effects to search for xallarap only and 
moreover, searching for a model including only xallarap effects does not provide significant improvements. For these reasons, 
we search for a solution that takes into account both orbital + terrestrial parallax and xallarap effects with a Markov Chain Monte 
Carlo algorithm (MCMC). We explore a large range of periods, from to 700 days, and find a.x^ improvement ix^ - 1717.7, 
Ax^ = -43) for periods above 250 days in comparison with the orbital parallax effects only. The;\f^ is essentially flat in the period 
range [250 - 500] days with a very shallow minimum around P - 290 days. 

The P = 290 days solution gives : = (0.958, -0.273), and thus a source orbital radius : - 6e - 3.74 AU. 

Kepler's third law (expressed in solar-system units). 



= P^M \M = M, + Mc 



and Newton's third law, 



imply 



M, 



Ml P2 [\+(MslMc)]^' 

From the position of the source relative to the red clump on the CMD diagram (Fig|2|l, we conclude that the source is a sub-giant 
situated in the bulge and, because the bulge is an old population, we infer that the source mass is close to a solar mass with 
an upper limit of 1 .2Mq. This mass limit and the long orbital period require a companion with > 7QMq, thus a black hole, 
which has an extremely low a priori probability. And if the companion is neither a black hole nor a neutron star, its mass has to 
be less or equal than the source mass since the source is an evolved star and a slightly heavier companion would therefore be 
much brighter. To explore these other possible star companions, we add a new constraint on the magnitude of the xallarap vector 
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1750 - 




150 200 250 

Pcriod(days) 

Fig. 5. a function of the period of the source'orbit from a MCMC algorithm with parallax and xallarap effects. The dashed 
line is the case without any constraint on the companion mass and leads to a black hole solution. The solid line is for a constraint 
< 2|^p2/3 ^ o.l8(P/yr)2/3) and leads to a solar mass companion with a minimum ofx^ at P = 170 days. 



in the MCMC program, assuming that M, < 1.2Mq and MdM^ < 1, which can be expressed as : 

3'/^ 

< ^^P'^' =0.18(^/3^)2/1 

The minimum of^^ (x^ = 1730) is obtained for a source orbital period equal to 170 days as shown in the Fig|5] When we put 
the corresponding parameters (P, 02, 62) in a differential-method program to reach a more accurate solution, we find;if2 - 1728.1. 

The xallarap vector of this solution {^e,N7 fo,£)=(-0.0142, 0.0940) implies a source orbital radius = Dg O^^e - OAQAU 
and a companion mass close to IMq. The MCMC algorithm permits us to explore an 1 1 -dimensional space (Iq, Ie, uq, P*, ^e,n, 
^E,E' ^E.N, ^E,E, P, ai, ^2)- We plot the Icr and 3cr limits of the [tteI = ^e as given in the Fig|6] The resulting parallax is then 
7Te = 0.94 ±6.10. 

3.4. Characteristics of the extended-source models with parallax and xallarap effects 

Considering the finite-source effects and parallax + xallarap effects, and the 16 observatories involved in the event monitoring, we 
have to fit 43 parameters (the 3 standard parameters, 1 for the angular size of the source, 2 for parallax, 5 for xallarap and 2x16 for 
the fluxes Fs and Ft of the different telescopes). The best ESPL fit model including parallax and xallarap effects (x^ = 1717.7) 
corresponds to a binary system in which the source companion is a black hole (see § 13. 3. 3b . One more reasonable solution could 
be the solar mass companion obtained when using a constraint on the xallarap (see § |3.3.3l l. This solution has x^ = 1728.1 for 
1745 data points and 43 fit parameters, to give x^/d.o.f - 1 .01, while the best ESPL fit with parallax effects only has;^'^ = 1760.5 
and the one without any parallax nor xallarap effects gives = 1995.4, a difference of Ax^ - 267.3. The corresponding best-fit 
parameters and their errors as determined from the light curve of three different models are shown in Table 1 and Table 2 (see 
also FiglT]). 
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ttE 

Fig. 6. as a function of the magnitude of the parallax vector he from MCMC runs including the constraint < 
(0.3 '''V3.7)p2/3 = 0.18(f'/yr)2/3 on the companion mass. ;r£ = 0.94 ± 0.10. 



TABLE 1 

OGLE-2007-BLG-050 Fit Parameters 
Extended-Source Point-Lens with Parallax 



Parameters 


Without 


Xallarap 


Xallarap 




xallarap 


Black hole 


Solar mass 






P = 290 days 


P = 170 days 


x' 


1760.5 


1717.7 


1728.1 


t(j {days) 


4221.9726 


4221.9725 


4221.9725 


cr,„ (days) 


0.0001 


0.0001 


0.0001 


Mo 


0.00204 


0.00215 


0.00214 


o-«o 


0.00002 


0.00003 


0.00003 


tE (days) 


68.09 


64.96 


65.11 


(Tt^ {days) 


0.66 


0.75 


0.75 


p. 


0.00450 


0.00473 


0.00471 




0.00004 


0.00006 


0.00006 


T^E 


0.12 


2.33 


0.94 




0.03 


0.07 


0.10 




/ 


1.00 


0.17 




/ 


0.06 


0.07 



Table 1. Fit parameters for three different models : 1/ with orbital parallax effects only, 2/ with orbital parallax -I- xaUarap 
(black-hole source companion), 3/ with orbital parallax -i- xallarap (solar-mass source companion). 
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TABLE 2 

OGLE-2007-BLG-050 Flux Parameters 
Extended-Source Point-Lens 
Orbital Parallax + Xallarap (solar mass companion) 



Observatory (Tf^ Ft, CTfi^ 



OGLE/ 


0.96 


0.01 


0.28 


0.01 


MOA/ 


0.96 


O.OI 


0.29 


0.01 


//FUN R New-Zealand (Auckland) 


1.08 


0.01 


3.06 


0.05 


//FUN R New-Zealand (Farm Cove) 


0.91 


0.009 


-27.0 


0.2 


//FUN / Ai-izona (MDM) 


0.99 


0.17 


242.3 


75.0 


//FUN R South Africa (Bronberg) 


0.232 


0.002 


0.28 


0.04 


//FUN / Chile (CTIO SMARTS) 


6.24 


0.06 


-3.9 


0.2 


//FUN / Arizona (Mt Lemmon) 


3.94 


0.04 


21.8 


0.1 


PLANET / South Africa (SAAO) 


5.60 


0.05 


-1.9 


1.6 


PLANET / Australia, Tasmania (UTas) 


2.74 


0.03 


163.8 


1.2 


PLANET / Chile (Danish) 


14.45 


0.15 


13.0 


0.4 


PLANET / Australia (Perth) 


0.793 


0.008 


6.5 


O.I 


PLANET / South Africa (Boyden) 


9.59 


0.09 


447.3 


0.6 


Robonet R Faulkes North (Hawaii) 


0.115 


0.001 


-0.26 


0.02 


Robonet R Faulkes South (Australia) 


0.129 


0.001 


-0.229 


0.009 


Robonet R Liveipool (Canaries Island) 


2.75 


0.03 


17.7 


0.3 



Table 2. Source flux and blending for telescopes fliat observed OGLE 07-BLG-050. The given values corresponds to the model 
with parallax and xallarap, in a case of a solar mass companion for the source. They do not change significantly for the other 
models. 



3.5. Lens mass and distance estimates 



Gouldl (119921) showed that if both 9e and tte could be measured, then the mass M and the lens-source relative parallax TTrei could 



be determined as given in Eq ll.2l and then the lens distance could be deduced from ; 

/ 1 1 

;r,-ei = lAU 

\Dl Ds 

The resulting characteristics of the lens are given in Table[3]for each model that we have presented : parallax only, parallax -H 
xallarap (black-hole companion) and parallax + xallarap (solar-mass companion). Due to the high parallax magnitude obtained 
with the "black hole" model (see Table [U, the lens mass is a brown dwarf (M = O.O25M0) in the extreme foreground (Dl - 
824/7c). Moreover, the extreme black hole mass (M, > TOM©), by itself, virtually rules out this model. We take this as evidence 
for unrecognized systematic errors at the Ax^ ~ 40 level, and hence do not beUeve inferences based on Ax^ at this level are 
robust. Systematic errors at this level are not uncommon for microlensing events. 

The model with a solar-mass companion is suspect as well, still with a brown-dwarf lens in the foreground, meaning that it 
results from the same systematics. We therefore conclude that the xallarap "signal" is probably spurious and we present these 
two models only for completeness. We expect that the presence of these systematics will corrupt the parallax measurements by 

of order ^ Av^ „ „ /Av^ „ ~ V(235 + 43/235) - 1 ~ 9%, which will impact the lens mass and relative parallax 

\ xaHarap+parallax' '^parallax i \ i ^ 7 x- x- 

estimates. However, this systematic error is too small to qualitatively impact the conclusions of this paper 

For the model with parallax eff'ects only, the lens star is a M-dwarf (Table O and situated in the disk, lying 5.5 kpc from the 
observer. With the added uncertainties due to systematics, the parallax becomes = 0.12 ± 0.03 ± 0.01, the lens mass estimates 
M = 0.50 ± 0. 14 Mq ( 28%) and the relative parallax nrei - 57.9 + 14.5yuas. For the rest of the analysis, we will only consider this 
model when the physical parameters of the lens are needed. 

As discussed bv . Ghosh et al. (2004) , future high-resolution astrometry could allow the direct measurement of the magnitude 
and direction of the lens-source relative proper motion fi and substantially reduce the parallax uncertainty and thus the stellar 
mass uncertainty. But according to our initial estimate of the relative proper motion (yu = 2.63 ± 0.08 mas/yr), it would take at 
least a 20 years to clearly detect the lens (especially since the source is very bright), but hopefully, within a decade, either ELT, 
GMT or TMT (giant telescopes) will be built, in which case the lens could be observed thereafter. 



4. Planet Detection efficiency 

4.1. Introduction and previous analyses 

To provide reliable abundance limits of Jupiter- to Earth-mass planets in our Galaxy, it is essential to evaluate the apparent 
non-planetary events, especially the well-covered high magnification events. A necessary step is to evaluate the confidence with 
which one can exclude potential planetary companions for each event. 
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TABLE 3 

OGLE-2007-BLG-050 Lens mass and distance 
Extended-Source Point-Lens with Parallax 



Parameters 


Without 


Xallarap 


Xallarap 




xallarap 


Black hole 


Solar mass 






P = 290 days 


P = 170 days 


Be 


0.48 ±0.01 


0.47 ± 0.041 


0.47 ±0.01 


M (MJ 


0.50 ±0.13 


0.025 ±0.001 


0.0618 ±0.0007 


TTrcl fytiasj 


58± 15 


1088 ±46 


440 ± 58 


Dl (kpc) 


5.47 ± 0.45 


0.82 ± 0.07 


1.77 ±0.20 



Table 3. Lens mass and distance for three different models : 1/ with orbital parallax effects only, 2/ with orbital parallax -H 
xallarap (black-hole source companion), 3/ with orbital parallax + xallarap (solar-mass source companion). 



Since OGLE-2007-BLG-050 presents strong finite-source effects, one may wonder whether a given planetary perturbation 
would have been so washed out by these effects as to become undetectable. Using many such efficiency calculations the aim is 
to determine the selection f unction to the underlying population of planets. 

Gaudi & Sackett ( 2000l) developed the first method to calculate detection efficiency for a single planet, which was extended 



to multiple planets detection efficiency by Gaudi et al. (2002), who analyzed 43 microlensing events from the 1995-1999 obser- 
vational seasons. Three of them were high magnification events [OGLE-1998-BLG-15 (A„„^ ~ 170), MACHO-1998-BLG-35 
(Amax ~ 100) and OGLE-1999-BLG-35 {A„,ax ~ 125)1. This 5-year a nalysis provided the first significant upper abundance limit 
of Jupiter- and Saturn-mass planets around M-dwarfs. iTsapras et al. 112003) and Snodgrass et al. (2004) derived constraints on 
Jovian planet abundance based on OGLE survey data of 1998-2000 and 2002 seasons respectively. 

Computing detection efficiency for individual events is thus required to estimate the frequency of planetary signatures in 
microlensing light curves, and a couple of complex eve nts have indeed be en analyzed separately. For example the high mag- 
nification event OGLE-2003-BLG-423 (A,„o., ~ 256) bv lYoo et alj (l2004bl) who found that the event was not as sensitive as it 
should have been if better monitored ov er the peak . Another high magnification {A,„a^ ~ 525) example is MOA-2003-BLG-32 / 
OGLE-2003-BLG-219 was analyzed bv lAbe et all (|2004) and Dong et al. (2006) (Appendix B). This well-covered event showed 
the best sensitivity to low-mass planets to date. Finall y, the highes t magn ification event ever analyzed, OGLE-2004-BLG-343, 
was unfortunately poorly monitored over its peak, and lPong et al. ( 20061) showed that it otherwise would have been extremely 
sensitive to low-mass planets. 



4.2. Planet detection efficiency in Einstein Units 

To characterize the planetary detection efficiency of OGLE-2007-BLG-050, we follow the lGaudi & Sackettl(l2000h method which 
consists of fitting binary models with the 3 binary parameters (d, q, a) held fixed and the single lens parameters allowed to vary. 
Here d is the planet- star separation in units o f Be, q the planet-lens mass ratio, and a the angle of the source trajectory relative 
to the binary axis. In lGaudi & Sackett ( 2000l) . the single lens parameters, uq, to and Ie, are related to a PSPL fit. In this analysis, 
we also fit the radius of the source p» (scaled to the Einstein radius) and compare the binary lens fits to the best ESPL fit for this 
event. 



From the resulting fitted binary lens xf^ ^ calculate the x^ improvement : Ax^^^^ ^ 



(d,q,a) 



' XeSPL' 



and Ayr, , is 

2 



compared with a threshold value I f Ax'^^ qa) ~Xc^ '■^^ ^' '^^ planetary (or binary) system is detected, while if Ax)^^ qa) " 
X^, it is excluded. iGaudi et al.l (12002!) argued that a threshold of 60 is high enough to be confident in excluding binary lens 
systems. 

For each (d, q), the fraction of angles < a < 27r that was excluded is called the "sensitivity" for that system. Indeed, the 
detection efficiency e(d, q) can be expressed as : 



e{d, q) 



- — 

2?r Jo 



@{Ax^{d,q,a)-xl)da 



where 9 is th e step funct i on. To perform the calculations of binary light curves, we use a binary-lens finite-source algorithm 
developed bv lDong et al. I (l2006h (Appendix A). The resulting grids of;^'^ as a function of d and a are shown in Fig |7] for some 
values of q. The complete computation has been done for every possible combinations between the following values of d, q and 
a : 



- q . 19 values with a constant logarithmic step over the range [10 ^, 10 ^]. 

- d:40 values with a constant logarithmic step over the range [0.1, 10]. 
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- a : 111 values linearly spaced from to 360°. 

The resulting detection efficiency diagram for OGLE-2007-BLG-050 is shown in FiglS] The first observation is that no 
planet is detected since there is no configuration that gives < -60. This event is very sensitive to the presence of planets, 
especially in the [0.8 - 1 .2] separation range in Einstein units, where the detection efficiency reaches 100% for Jupiter mass ratios 
{q - 9x10""*), 75% for Neptune mass ratios {q - 5x10"^) and 10% for Earth mass ratios (3x10^^). In larger separation ranges, as 
[0.4 - 2.7] R^, we exclude Jupiter mass ratios with 95% confidence. 

In future statistical analyses of microlensing planetary detection efficiency, one will likely be forced to use a higher exclusion 
threshold than 60 because, while planets can sometimes be reliably excluded at this threshold (as in the present case), it is 
unlikely that they can be reliably detected at this level, particularly in high-magnification events. Because we cannot predict the 
exact threshold that will be adopted by future studies, we show both our exclusion level (Ax^ > 60) and a somewhat arbitrarily 
chosen value, Ax^ > 250. The important point is that the detection efficiency diagrams in the two cases (Fig|8]and with a threshold 
equal to 250 in Fig|9]l are very similar. 



4.3. Planet detection efficiency in pliysical units 

Having an estimate of the angular Einstein radius 6e, the distance Di of the lens from the observer and the lens mass M, we 
derive estimates of the physical parameters ir±,mp) for the tested planetary models, where r± is the projected separation between 
the planet and its host star and nip the planet mass, and calculate the associated detection efficiency. 

r_L (AU) = (kpc) 0E (mas) 

OTp — qM 

To simplify the translation between efficiency diagrams in Einstein units and physical units, we have considered the values 
of M, Di and flg as perfectly known. A proper analysis would involve a convolution of the detection efficiency map in terms of 



native parameters d,q over the probability density distribution of the primary lens parameters (e.g. lYoo et al. I ( l2004bl) ). While 
our procedure of keeping M, Di and 6^ fixed is an approximation, considering the logarithmic scale of the efficiency maps, the 
uncertainties in the primary lens parameters will not have an important effect on the shape of the resulting efficiency diagrams. 

We take the parameters related to the fit with extended source and parallax effects, where M - 0.50 ± 0.14Mq, Di - 
5.47 + 0.45 kpc and 0^ - 0.48 + 0.01 mas. The resulting detection efficiency diagram in physical units is shown in Figl8]as well, 
but the corresponding axes are those on the top and the right of the graphic. This demonstrates that OGLE-2007-BLG-050 is 
sensitive to Neptune-mass planets as well as some Earth-mass configurations. Indeed, for a [1.8 - 3.1] AU projected separation 
range between the planet and the lens star, Jupiter, Neptune and Earth-like planets are excluded with a 100%, 95% and 10% 
confidence respectively. For a range of [1.4 - 4] AU, the detection efficiency reaches 100% for Jupiter mass planets and 75% for 
Neptune mass planets, and for a much bigger range of [0.6 - 10] AU, Jupiter-like planets are excluded with a 75% confidence. 



4.4. Planet detection efficiency as a function of central caustic size 



Chung et al.l (l2005h analyzed the properties of central caustics in planetary microlensing events in order to estimate the pertur- 
bation that they induce. They gave an expression for the central-caustic size as a function of the planet-star separation and the 
planet/star mass ratio. Several authors have considered the size and shape of the central caustic as a func tion of the parameters 
of the planet for high-magnification events (iGriest & Safizadeh"l998l IPominik |[l999t IPong et alJEo09bl) . In the analysis of the 
cool Jovian-mass planet MOA-2007-BLG-400Lb, [Bong et al. (2009b) conducted the initial parameter space search over a grid 
of {w,q) rather than id,q) where w is the "width" of the central caustic. For MOA-2007-BLG-400, the angular size of the central 
caustic is smaller t han that of the source (w/p ~ 0.4), and w can be directly estimated by inspecting the light curve features. 
Dong et al. I (l2009bh find the (w q) parametrization is more regularly defined and more efficient in searching parameter space than 
(d,q). 

The source size of OGLE-2007-BLG-050 is p = 0.0045 which is relatively big, and since finite-source effects smear out 
the sharp magnification pattern produced by the central caustics, one way to present the planetary detection efficiency results 
is to estimate the ratio w/p that is reached at the detection/exclusion limits. Assuming that detectable planets should produce 
signals > 5%. lHan & KimI (12009 ) estimated the ratio w/p must be at least equal to 0.25. Here we present the detection efficiency 
diagram in (d,w/p) space in FigHOl still considering Ax^ > 60 as the criterion of exclusion. This diagram shows a clear frontier 
in red at w/p values between 0.1 and 0.3 above which the detection efficiency is greater than 50%, which also corresponds to 
the 50% detection's contours in Fig[8] On this frontier, the value of w/p goes down to 0.1 for d ~ I and increas es to 0.3 for 
d » 1 or d « 1. Our realistic estimate of detection efficiency is in general agreement with the simple criterion in Han & KimI 
(l2009h. ' jiven the hig h photometric precision and dense sampling, our data allow detections below the 5% threshold adopted 
bv lHan&KimI ( |2009|) . We also note that the w/p threshold is weakly dependant on d, which is a result of the enhancement in 



detection efficiency of the resonant caustics at small mass ratios. 
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We have presented a new way of visualizing the detection efficiency in (d,w/p) space. It offers a physically straightforward 
way to understand the planetary sensitivity in events with pronounced finite-source effects. We find that the data obtained by 
current observation campaigns can probe planetary central caustics as small as ~ 20% of the source size for high-magnification 
events. 

5. Conclusion 

OGLE-2007-BLG-050 is a rare case of a high magnification event with well measured finite source effects and detectable parallax 
effects. This leads to an estimate of the angular Einstein radius Oe = 0.48 ± 0.01 mas, the parallax tte = 0.12 ± 0.03, the mass 
M = 0.50 ± 0.14 Mq and distance Dl = 5.5 ± 0.4 kpc of the lens star. This is only the second reasonably precise mass estimate 
(to within 28%) for an unseen single object using any method. 

When computing planet detection efficiency, we did not find any planetary signature and the resulting maps in (d, q, a), where 
d is the planet-star separation in Einstein units, q the planet-lens mass ratio, and a the angle of the source trajectory relative to the 
binary axis, reveal a good sensitivity to low mass ratios q, with a 75% and 10% efficiencies for Neptune- and Earth-mass ratios 
respectively in the range [0.8 - 1.2] R^,, and a 100% detection efficiency for Jupiter-mass ratio in [0.4 - 2.7] Re- 

It also permits the calculation of efficiency maps in physical space ir±,mp), where r± is the projected planet/star separation 
and nip is the planet mass. Here we show that this microlensing event is very sensitive to Neptune-mass planets and has (10%) 
sensitivity to Earth-mass planets within a [1.8 - 3.1] AU projected separation range. 
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Fig. 7. Binary-lens finite-source grids of as a function of (x, y) where x = d cos a and y = dsina for different values of q. The 
value appearing in the upper part of each diagram corresponds to the value of log q. The color scale shows the variations of the 
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Fig. 8. Resulting detection efficiency diagram for d and q ranges of [0.1 - 10]/?e and [10^^ - 10"^]M;f.„.s and detection efficiency 
diagram in physical units (r±,mp) if considering the upper and right axes. The contours indicate 25%, 50%, 75% and 90% 
efficiency, with an excluding threshold equal to 60. 
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Fig. 9. Same as Fig[8] except with a hypothetical threshold of Ax^ - 250. Comparison with Fig|8] shows that planet sensitivity 
does not depend strongly on threshold. 
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Fig. 10. Resulting detection efficiency diagram in {d,w/p) space. This diagram shows a clear frontier in red atw/p values between 
0.1 and 0.3 above which the detection efficiency is greater than 50%. This frontier corresponds to the 50% detection's contours 
in the Figi] 
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